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1.  INTRODUCTION 


Injury  to  the  gastro-intestinal  tract  is  one  of  the  primary  biological  effects  of  strong  blast 

>  exposure.  The  threshold  for  injury  appears  to  lie  between  that  for  the  larynx  and  the  lung  and 
therefore  may  be  a  critical  element  in  defining  Damage  Risk  Criteria  for  occupational  exposure. 
In  combat  circumstances,  this  injury  may  be  debilitating. 

Little  is  known  about  the  mechanics  of  the  injury  process  or  its  relation  to  external  blast 

>  characteristics.  In  fact,  little  is  known  about  the  pressure  environment  within  the  abdomen  during 
blast  loading.  Laboratory  experiments  at  JAYCOR  using  perfused,  small  animal  intestine  have 
provided  most  of  our  current  understanding  of  the  phenomena  [1]. 

In  those  tests,  the  static  bursting  strength  of  various  sections  of  rabbit  intestine  were 
determined.  There  seems  to  be  little  systematic  variation  with  location  along  the  gastro-intestinal 
*  tract  or  with  animal  weight,  age,  or  sex.  The  material  properties  were  in  agreement  with  Yamada 

[2].  The  most  prevalent  form  of  injury  is  hemorrhage  without  associated  rupture,  so  that  it  is  likely 
that  threshold  injury  is  controlled  by  the  failure  of  the  vascular  system. 

Dynamic  experiments  were  then  conducted  in  which  isolated  sections  of  intestine  were 

>  placed  in  a  water  filled  chamber  and  exposed  to  an  external  blast  loading.  Damage  to  the  intestine 
could  be  directly  correlated  with  the  location  of  air  bubbles  within  the  section.  Violent  motion  of 
the  tissue  was  observed  and,  during  the  overexpansion  phase,  rupture  sometimes  occurred. 

Tests  made  with  balloons  filled  with  air  and  water  and  immersed  in  water  produced  a  similar 
dynamic  response.  High  speed  movies  showed  surface  deformations  that  were  primarily  limited  to 
the  local  region  at  the  bubble.  As  the  viscosity  of  the  interior  fluid  was  increased,  the  deformations 
became  even  more  localized.  Pressure  probes  placed  inside  the  bubble  showed  that  the  surface 
motion  could  be  correlated  with  the  variations  in  the  pressure  time  histories.  This  correlation 
suggested  that  bubble  pressure  measurements  could  be  used  to  understand  the  injury  process. 

>  Interpretation  of  the  observed  damage  in  isolated  intestine  sections  was  clouded  by  the 
differences  from  the  in  vivo  situation.  The  lack  of  a  blood  supply  causes  the  intestine  to 
deteriorate  rapidly  and,  of  course,  eliminates  the  possibility  of  observing  hemorrhage. 

A  surgical  technique  was  developed  by  Yu  and  Vasel  that  allowed  an  isolated  rabbit 
(  intestine  to  be  oxygenated  by  the  animal’s  own  cardiopulmonary  system.  The  perfusion  results 

from  connecting  the  abdominal  aorta  to  the  cranial  mesenteric  artery  and  connecting  the  portal 
vein  to  the  caudal  vena  cava.  The  intestine  can  be  maintained  in  a  nearly  in  vivo  state  while  being 
subjected  to  blast  loading  in  a  water  chamber. 
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Results  obtained  with  the  previous  isolated  preparation  were  confirmed  and  thresholds  of 
hemorrhagic  injury  at  lower  levels  were  observed.  Occurrence  of  injury  was  not  completely 
correlated  with  the  peak  chamber  pressure,  however,  suggesting  that  local  conditions  were 
important. 

Direct  measurement  of  the  stress  in  the  wall  section  proved  to  interfere  with  the  dynamical 
motion  and  a  noninvasive  technique  was  sought.  Pressure  transducers  placed  within  the  bubble 
showed  the  same  behavior  as  had  been  observed  in  the  isolated  preparations.  Moreover,  the 
pressure  differential  between  the  bubble  and  a  position  in  the  fluid  near  the  bubble  seemed  to 
correlate  better  with  injury. 

As  a  result  of  these  findings,  a  relation  between  air  bubbles  and  injury  was  proposed.  The 
pressure  differential  across  the  membrane  gives  a  measure  of  the  dynamics,  hence  the  stress  in  the 
tissue.  An  analogy  was  drawn  with  the  relation  between  surface  tension  and  pressure  difference  in 
a  bubble.  The  amount  of  test  data  on  which  this  concept  is  based  is  small,  however,  and  there  are 
many  uncontrolled  factors  in  dealing  with  the  entire  gastro-intestinal  tract.  A  simpler  circum¬ 
stance  was  sought  in  which  to  develop  and  test  these  ideas. 


2.  SURROGATE  EXPERIMENTS 


The  objective  of  the  experiments  is  to  provide  data  on  gas  and  lumen  dynamics  under 
controlled  circumstances  so  that  the  hypotheses  proposed  can  be  tested.  At  the  same  time,  we 
want  to  perfect  techniques  for  measuring  the  tissue  dynamics  that  can  be  used  in  future  animal 
tests. 


Lamb  cecum  sections,  in  a  commercially  available  and  preserved  form,  were  used  as  the 
surrogate  material.  Tests  were  run  to  demonstrate  that  they  had  material  strengths  comparable  to 
that  of  rabbit  intestine.  A  fixture  was  designed  to  be  placed  in  the  water  chamber,  holding  the 
section  at  a  circular  ring.  The  fixture  formed  an  enclosure  surrounding  the  membrane  with  flow 
holes  at  the  side.  See  Figure  1.  By  adjusting  the  volume  of  water  and  gas  under  the  membrane, 
the  curvature  could  be  varied  continuously  from  convex  to  concave  and  the  bubble  could  be  given 
any  volume. 

Pressure  probes  were  mounted  at  fixed  locations  on  either  side  of  the  initial  membrane 
position.  Pressure  measurements  were  also  taken  in  the  chamber,  away  from  the  fixture. 

As  the  intensity  of  the  blast  was  increased,  as  measured  by  the  peak  chamber  pressure,  the 
maximum  pressure  reached  inside  the  bubble  and  the  frequency  of  bubble  oscillation  increased. 

The  differential  pressure  between  the  inside  and  outside  probes  also  increased  with 
intensity.  For  positively  curved  membranes  (bubble  on  the  concave  side)  the  differential  pressure 
first  goes  negative  before  rising  to  its  positive  maximum.  For  negative  curvature  (bubble  on  the 
convex  side),  the  curve  reverses  in  sign.  Rupture  of  the  membrane  usually  occurred  in  this  case. 

The  differential  pressure  maximum  did  not  vary  significantly  with  the  size  of  the  bubble,  but 
it  decreased  as  the  membrane  was  stretched  flatter.  The  differential  pressure  also  decreased  as  the 
initial  bubble  pressure  was  increased. 

These  observations  form  a  basis  on  which  to  develop  and  test  the  analytical  model  discussed 
in  tf  e  next  section. 


4  MOUNTING  BOLTS  TRANS  A. 


STAND-OFF  NUTS  FOR 
TRANS.  CABLE  PROTECTION 


Figure  1 .  Schematic  of  single  curvature  membrane  test  fixture. 


ANALYTICAL  MODEL 


DYNAMICS  OF  THE  LIQUID  MOTION 

We  use  the  statement  of  energy  conservation  to  develop  a  dynamical  equation  of  motion 
The  rate  of  change  of  kinetic  energy  in  the  liquid  outside  the  membrane,  E,  is  equal  to  the  work 
done  by  pressure  forces  at  the  interface  and  at  the  far  boundares  less  the  rate  at  which  energy  is 
dissipated  within  the  liquid  flow. 


dE 

dt 


(pAv).  -  (pAv)  -  aE 
1  D 


(i) 


where  the  kinetic  energy  is  given  by 


...  1  2 
E  -  JJT  2  PL  V  dV 


(2) 


A  dissipation  rate  proportional  to  the  kinetic  energy,  as  occurs  in  turbulent  flow,  has  I'een 
assumed.  The  pressure  at  the  interface  on  the  liquid  side  is  equal  to  the  pressure  in  the  gas,  pq, 
less  the  differential  pressure  caused  by  the  tension  in  membrane,  6pM, 


Pi  =  PG  ”  5PM 


(3) 


The  contribution  of  the  membrane  depends  on  its  material  properties  and  its  deformation. 

We  imagine  a  coordinate  system  based  on  the  liquid  in  which  one  coordinate  direction,  x,  is 
parallel  everywhere  to  the  flow.  The  flow  area  normal  to  the  coordinate  x  is  designated  A(x).  The 
liquid  velocity  is  then  always  normal  to  A  so  that  the  statement  of  mass  conservation  becomes 

PL  v(x)  A(x)  =  constant  .  (4) 

Since  the  liquid  is  assumed  to  be  incompressible  and  since  the  geometries  must  be 
conformal  to  the  membrane  position,  then 

v(jc)  A(x)  =  vjAj  (5) 

The  differential  volume  in  this  coordinate  system  is  given  by 

dV  =  A  dx  .  (6) 


Combining  the  above  relations  gives 
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The  rate  of  change  of  the  energy  is  then 
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The  rate  of  change  of  the  volume  flux  is 


dv.  dA 

—  (v  A.)  -  A,  — -  +  v.  - — 
dt  i  i  i  dt  l  dt 


dv.  f dx. 
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Combining  (7)-(9)  with  (1)  yields  the  equation  of  motion 


dv 

1  dT  +  R  Vi  “  (PG  ~  5pM  “  PB)/PL  ■  QlVi 


(10) 
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where  pg  is  the  applied  pressure  at  the  far  boundary  and  the  inertia  and  kinematic  coefficients  are 
given  by 


I  -  A 


[XB  dx 

i  J  A(x) 


fdAj  I^B  dx  _  1 
[dxj .  J  A(x)  2 


I  X. 
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DYNAMICS  OF  THE  BUBBLE  MOTION 

The  rate  of  change  of  the  bubble  volume  is  equal  to  the  volume  flux  in  the  liquid 


— - A.v. 

dt  li 


since  the  liquid  is  incompressible.  Because  the  blast  loadings  are  so  rapid,  we  assume  that  process 
is  adiabatic  within  the  bubble  so  that  for  an  ideal  gas 

PGVb  “  PG(0>  Vb(0)  '  (14) 

Finally,  the  pressure  differential  due  to  the  membrane  must  be  found  from  experiment.  For  an 
elastic  material 


%  '  <PM<V 


GEOMETRIC  PARAMETERIZATION 

Various  flow  geometries  can  be  simulated  through  the  choice  of  V(x).  For  example,  a 
spherical  bubble  in  an  infinite  medium  would  be  characterized  by 

V(x)  -  |  n  x3  (16) 


XB-“  * 
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Using  (6),  (11),  and  (12)  yields 


A(x)  -  4xx 
I  -  x 


R-| 


which  produces  Rayleigh’s  equation  [3] 


dv  3  2  G  ® 

x  —  +  —  v  -  - 

dt  2  p. 


Similarly,  a  plug  of  water  of  radius  r  and  length  L  yields 


V(x)  -  jrr  x 


xfi  -  x.  +  r 


so  that 


A(x)  -  nr 
I  -  L 


and  the  dynamical  equation 


T  dVi  12  Pa  “  PL 

L  5T  -  2  vi  “  — TT 


(18) 


(19) 

(20) 


(21) 
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For  the  surrogate  experiments  described  earlier,  the  appropriate  geometry  is  a  spherical 
sector  anchored  at  a  ring  of  radius  rj  with  an  outer  boundary  at  a  radius  X2 ■  Then 


which  yields 


V (x)  -  |  x^3r2  +  x2) 
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LIQUID  PRESSURE  FIELD 

Because  the  surrogate  experiments  collected  pressure  data  at  points  within  the  liquid  it  is 
desirable  to  form  an  expression  for  the  pressure  field.  If  the  previous  derivation  were  carried  out 
for  a  region  of  liquid  bound  by  x  and  Xjj,  then  we  would  have  found 


rB  dx  d_ 
A(x)  dt 


'V*)  '  <p  ■  V  +  \ 


Using  (5)  and  rearranging  the  left-hand  side  of  (24)  yields 


1(x)  Tt  ('lvA)  -  -  pb)  +  I  v 
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while  the  same  form  at  the  interface  is 


ir  h  ‘-ViV  -  (pg  -  5pm  -  V  + 1 

1 


Taking  the  ratio  of  (25)  and  (26)  and  solving  for  p  yields 


A. I (x) 


.  ,  1  2  f  1  2}  i  v 
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(26) 


(27) 


NUMERICAL  SOLUTION  ALGORITHM 

The  set  of  coupled  ordinary  differential  equations  for  the  position  and  velocity  of  the 
membrane  coordinate  is  solved  by  a  fifth-order  Runge-Kutta  algorithm  [4].  Time  steps  were  taken 
to  be  0.05  msec  and  variations  were  made  to  assure  that  accurate  solutions  were  being  obtained. 


4.  MODEL  VALIDATION 


PARAMETER  ASSIGNMENTS 

For  the  surrogate  test  fixture  described  earlier,  the  following  geometric  parameters  apply. 
The  5  cc  bubble  forms  a  spherical  section  with  base  radius  of  1.68  cm  and  height  1.00  cm.  Based 
on  the  previous  experimental  observations,  we  assume  that  the  membrane  distorts  only  in  the 
region  of  the  bubble  and  maintains  the  geometry  of  a  spherical  section.  The  partly  enclosed  fixture 
has  an  effective  radius  of  1.91  cm  as  measured  from  the  bubble  center.  We  then  can  use  the 
geometry  described  by  (22)  and  (23)  with  the  values 

rj  =  1.68  cm 

T2  =  1.91  cm 

x(0)  =  1.00  cm 

The  differential  pressure  in  the  membrane  is  taken  from  tests  data  to  increase  at  a  rate  of  3 
psi/cm.  The  flow  is  assumed  to  have  no  damping. 

For  the  negative  curvature  runs,  we  have 
r  i  =  1.91  cm 
T2  =  1.91  cm 
x(0)  =  -1.91  cm 

The  chamber  pressure  is  taken  from  the  measured  trace  for  each  shot  and  is  characterized 
by  a  peak  value  and  duration  to  the  first  zero  crossing.  The  durations  varied  from  4.4  to  5.0  ms. 


DISCUSSION  OF  RESULTS 

Figure  2  compares  the  bubble  pressure  time  histories  for  a  series  of  tests  in  which  the 
chamber  pressured  varied  from  27  to  116  psig.  The  plots  are  not  to  scale  so  only  shape  can  be 
compared  here.  There  seems  to  be  generally  good  agreement  on  the  number  and  shape  of  peaks. 
The  increasing  oscillation  frequency  with  blast  intensity  is  captured. 

Figure  3  shows  the  quantitative  comparison  of  these  results.  The  decrease  with  intensity  of 
the  time  between  peaks  is  accurately  predicted  by  the  model  as  is  the  first  maximum  and  minimum 
in  the  trace.  At  high  intensities,  the  calculation  tends  to  over  predict  the  second  and  third  maxima. 
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Figure  3.  Comparison  of  calculation  and  measurement  of  bubble  pressure  extrema 
parameters  for  the  conditions  shown  in  Figure  2. 
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A  review  of  the  experimental  data  revealed  that  the  membrane  struck  the  inside  transducer  in 
these  cases.  The  impact  may  have  slowed  the  collapse  process  and  caused  a  smaller  peak  to  be 
reached. 

Figure  4  compares  the  bubble  pressure  time  histories  for  the  cases  of  positive  and  negative 
curvature  in  the  membrane.  The  change  in  frequency  is  captured.  Figure  5  shows  the  calculated 
pressure  differential  between  the  inside  and  outside  transducers.  The  dramatic  reversal  of  the 
trace  reported  in  the  experimental  data  is  seen.  Figure  6  presents  the  comparison  of  the 
calculated  and  measured  pressure  differential  maximum  values.  Agreement  is  only  qualitative  in 
this  case.  Again,  the  impact  with  the  inner  transducer  may  have  affected  the  results.  However, 
both  5  and  10  cc  bubbles  produced  nearly  the  same  differential  pressures,  confirming  the 
experimental  observations. 

The  predicted  bubble  responses  did  not  vary  significantly  with  the  assumed  stiffness  of  the 
membrane  for  the  range  of  values  measured,  indicating  that  the  pressure  differentials  are 
primarily  controlled  by  the  flow  dynamics  in  the  water.  Another  evidence  of  this  dependence  was 
that  the  response,  especially  the  pressure  differential,  varied  with  the  assumed  flow  geometry. 
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Figure  4.  Comparison  of  measured  and  calculated  bubble  pressure  time  histories  for 
cases  when  the  membrane  has  positive  and  negative  curvatures. 
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Figure  5.  Comparison  of  calculated  differential  pressure  time  histories  for  positive  and 
negative  curvature  of  the  membrane. 
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Figure  6.  Comparison  of  calculated  and  measured  differential  pressure  maximum 
values. 
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5.  CONCLUSIONS 


The  analytical  model  provides  a  reasonably  good  explanation  of  the  observed  dynamics  of 
the  bubble.  In  future  tests,  it  will  be  used  to  predict  the  response  in  a  complex  wave  environment 
and  to  anticipate  membrane  rupture.  Eventually,  these  concepts  will  be  employed  in  developing  an 
understanding  of  the  mechanisms  of  gastro-intestinal  injury  in  vivo. 
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